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Abstract. 

We present results from a fully relativistic investigation of matter fallback in a 
supernova, assuming that the central compact star is a black hole (BH). Our goal is 
to answer the following questions: can we infer from the light curve whether a stellar 
black hole (BH) has formed in the aftermath of a supernova explosion ? Did SN1897A 
produce a BH and, if so, when will the hole become detectable ? We show that it 
should not be possible to infer the presence of a BH in the remnant for another several 
hundred to a thousand years. In the case of SN1987A, we estimate the BH luminosity 
to be L ~ 5 x 10 34 erg s _1 , which is well below the present day bolometric luminosity 
of the remnant (~ 10 36 ergs -1 ; Suntzeff 1997 [1]). 

INTRODUCTION 

Stellar evolution calculations show that stars with a main sequence mass in the 
range 8-19 M Q finish their lives with a compact core of ~ 1.4M (Woosley & 
Weaver 1995 [2]). For these stars the amount of material that falls back toward 
the core in the aftermath of its collapse is negligible, so it is likely that they give 
birth to neutron stars. The fate of stars with main sequence mass in the range 
19-25 M & is far less obvious. At the end of their evolution, they have core masses 
around 1.6-1.8 M Q and a variable amount of matter, 0.1-0.3 M & may fall back. 
Typically, these explosions leave a central compact remnant of 1.7-2.1 M© whose 
nature depends critically on the maximum (gravitational) mass of a neutron star, 
M crit , and hence on the equation of state of dense nuclear matter. To date, the 
value of M crit is quite uncertain: 1.5M Q < M crit < 2.5M for a nonrotating star. 

The mass of the progenitor of SN1987A (M = 18-21 M Q ) falls in the range of 
masses where the outcome of core collapse is uncertain; from a theoretical point of 
view, we do not know if a neutron star or a black hole formed during the explosion. 
However, if a compact object of any kind is present, it is probably accreting from 



the progenitor stellar material (Chevalier 1989 [3]). The bolometric light curve 
observed to date can be explained by the "standard" theory of Type II supernovae. 
It is therefore important to determine what extra luminosity an accreting central 
component would produce and when its presence might be discernible. In this 
paper we present results from the first fully relativistic investigation of supernova 
fallback in presence of a BH including the effects of hydrodynamics and radiation 
transport. 

METHOD AND NUMERICAL RESULTS 

We constructed a spherically symmetric, general-relativistic radiation hydrody- 
namic Lagrangian code capable of handling the transfer of radiation from the early 
phase when photons diffuse through a high temperature, expanding cloud, to the 
late stage when the hydrogen envelope has recombined and most of the ejecta are 
transparent (Zampieri et al. (1997) [4]). We consider fallback from a hydrogen 
"cloud" ; all of the emitted energy comes from the release of heat residing in the gas 
originally, or is generated by compressional heating in the course of accretion onto 
the BH. Radioactive energy sources have not been included in our calculations. At 
the onset of evolution the cloud has homogenous density po and a velocity profile 
which scales linearly with radius r. The temperature profile is taken to be equal to 
the "radiative zero solution" given by Arnett (1980) [5] and the degree of ionization 
is computed throughout evolution according to the Saha equation. Four parameters 
uniquely specify the state of the "cloud" : the total mass M c i ou d, the radius r out , the 
sound speed at the inner boundary c S; o, and the ratio k of the accretion timescale 
ta,o = GM bh /c^ (M bh mass of the BH) to the expansion time t = r out /V (V 
maximum expansion velocity of the ejecta). 

We find two phases of evolution. In the first stage, gas expands adiabatically 
in the outer region and accretes close to the central hole. During this phase, the 
cloud is radiation dominated. A second phase follows, during which hydrogen 
recombines and the "cloud" becomes optically thin. Fig. 1 shows the temperature 
profile of model A characterized by M doud = 1M (p = 2 x 1CT 5 g cm -3 ), r out = 
3 x 10 12 cm, k = 0.1 and c Sj0 = 10 8 cm s -1 . A recombination wave develops 
that propagates rapidly through the envelope. The recombination front stalls as 
soon as it approaches the innermost region where compressional heating due to 
accretion overcomes radiative cooling. Outside the recombination front the flow is 
transparent but it remains optically thick inside for a long time. The density 
and velocity profiles display a self-similar character and are unaffected by the 
propagation of the recombination front. 

As recombination begins, the huge amount of internal energy released gives rise 
to a big bump in the light curve at t ~ 10 days, as illustrated in Fig. 2. After the 
emission maximum, the luminosity falls off abruptly. During the late evolution, the 
light curve is entirely powered by accretion onto the central BH and the luminosity 
decreases as a power law with time. At late times, the evolution proceeds through 
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FIGURE 1. Gas temperature T versus radius r (in units of the initial accretion radius 
r a ,o = GMbh/c 2 s .o) a t selected times for model A. From the top to the bottom: t/t a _ = 0, 
729, 3757, 4755, 4870, 11175. 

a sequence of quasistationary states, and m (in Eddington units) can be estimated 
using the scaling relation for a pressureless fluid (Colpi et al. 1996 [6]; equation 
[29]) 

m(t) ~ [47T 2 / 3 /9]p t k es c (t/t )~ 5/3 , (1) 

where k es is the electron scattering opacity. The luminosity / (in Eddington units) 
is close to that derived by Blondin (1986) [7] for stationary accretion 

I~3x lCT 7 (M bh /M & y 1/3 m 5/6 . (2) 

In Fig. 2, the dotted line at the bottom denotes the late-time accretion luminos- 
ity estimated using equations (1) and (2): the analytic extrapolation is in good 
agreement with the computed light curve. 

THE EMERGENCE OF A BH IN SN1987A 

Slightly after maximum, SN1987A was powered by the radioactive decay of 56 Co 
and 57 Co, and at present its light curve is consistent with emission from the decay 
of 44 Ti (Suntzeff 1997 [1]; Fig. 2). Can we discern the luminosity emitted by the 
accreting BH above the contribution resulting from 44 Ti ? An upper limit to the 
accretion luminosity of SN1987A can be inferred from equations (1) and (2) (for 
details see Zampieri et al. (1997) [4]). Following Chevalier (1989) [3], we adopt 
p ~ 10 9 cgs and t ~ 7000 seconds, yielding a luminosity 




8 x 1(T 3 
~ (M 6 ,/M ) 1 /3[ t (years)] 25 / 18 

Equation (3) implies L ~ 5 x 10 34 ergs _1 after t ~ 10 years, which is well below the 
present day bolometric luminosity of the remnant (~ 10 36 ergs _1 ; Suntzeff 1997 [1]) 
and also smaller than the luminosity estimated to result from radioactive decay. 
Thus, there is no observation that rules out the possibility that a BH resides inside 
the SN1987A remnant. In Fig. 2 we plot the late-time light curve of SN1987A 
calculated using (3). As radioactive decay plummets at around 270-2700 years, it 
is clear that the BH would only appear after about 900 years irrespective of the 
detailed numbers. After this time has elapsed, the luminosity of the remnant would 
be ~ 10 32 ergs _1 , too dim for detection with present technology. The inner part of 
the accretion flow will be still optically thick and will be emitting roughly a black 
body spectrum peaked in the visible band. 
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FIGURE 2. The Circles give the solution for model A. The Dotted line is the extrapolation of 
the late-time evolution. The Triangles are the bolometric luminosity of SN1987A. The Dashed 
lines represent the contribution from the decay of radioactive elements (0.07 M@ of 56 Co and 
~ 5 x 10 _5 M Q of 44 Ti). Finally the upper dotted line denotes the expected bolometric luminosity 
emitted by a putative BH in SN1987A (equation [3]). The arrow marks the time (t ~ 900 years) 
of the BH emergence. 



